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1 .  INTRODUCTION 


The  electromagnetic  pulse  (EMP)  analysis  of  a  system  frequently  focuses  on 
the  response  of  cables  excited  by  the  transient  fields.  A  large  variety  of 
analysis  tools  and  techniques  is  available  to  aid  in  the  study  of  the  cable 
response.  EMP  simulators  are  available  for  experimental  studies  of  system 
response.  Difficulties  arise,  however,  when  we  attempt  to  incorporate  the 
results  of  cable  coupling  codes  into  circuit  analysis  codes  for  damage  assess¬ 
ments.  Similarly,  obtaining  threat-level  responses  with  typical  simulators  is 
possible  for  only  very  short  cables  and  compact  systems.  This  report  presents 
a  solution  to  these  problems  using  an  equivalent  source  representation  of  the 
cable.  Initially,  the  technique  is  developed  for  a  single  coaxially  shielded 
cable.  Then  the  approach  is  extended  to  multi  conductor  cables. 

2.  RESPONSE  OF  A  COAXIALLY  SHIELDED  CABLE 

The  response  of  a  coaxial  cable  due  to  a  distributed  excitation  is  re¬ 
viewed  in  section  2.1.  An  equivalent  and  more  compact  source  representation 
is  developed  in  section  2.2.  Section  2.3  investigates  the  validity  of  the 
equivalent  representation  when  the  line  is  loaded  with  nonlinear  terminations. 

2. 1  Distributed  Excitation 

The  response  of  a  coaxially  shielded  cable  with  distributed  excita¬ 
tion  may  be  treated  like  the  two-wire  line  problem  shown  in  figure  1  .  The 
distributed  voltage  source,  Vs(x,ai),  is  related  to  the  complex  transfer  imped¬ 
ance  of  the  shield  and  the  external  current  flowing  on  it.  The  distributed 
current  source,  Ia(xfw),  accounts  for  coupling  through  the  holes  of  the 
shield.  Detailed  explanations  of  the  shield  coupling  mechanisms  and  solutions 
of  the  coupling  problem  are  provided  in  numerous  reports  (e.g.,  [1],  [2], 
[3D*  For  this  paper  it  is  sufficient  to  define  the  general  coupling  problem 
as  two  arbitrary  distributed  sources — a  series  voltage  source  and  a  shunt 
current  source.  The  transmission  line  of  figure  1  has  a  characteristic  imped¬ 
ance,  K,  a  propagation  characteristic,  r,  and  linear  terminations  Z 1  and  Z2  at 
x  -  0,1,  respectively. 

Schelkunoff  [4]  provides  a  solution  for  the  terminal  response  of 
the  transmission  line  shown  in  figure  1 : 

1  1 

V(0)  =  Jv  (5)Gvl(0£)d5  +  J  Is(S)Gv2(0£)dS 

o  o 

K Z.  r 

=  -  J  V.®  {K  cosh  TO  -  $  +  Z2  sinh  T  (1  -  $)}d£  - 


•References  are  listed  at  end  of  report. 


K  Z  f 

-  J  Is  (%)  {K  sinh  TO  -  %)  +  Z2  cosh  T (1  -  $)}d$ 

O 

1  1 

V(1)  =  J  v  (5)GV1  (1,1) di  +  Ji,  ©  gV2  0,5)d5 

o  O 

KZ,  J 

=  I  V(^)  {K  cosh  +  Zj  sinh  T^}d^ 

o 

k2z  1 

-  -^1  j  Is  (4)  {K  sinh  +  Zj  cosh  r^}d^ 

o 

where 


Figure  1.  Transmission  line  with  distributed  excitation. 


The  line  response  may  be  computed  using  the  proper  frequency  de¬ 
scriptions  of  the  line  parameters,  terminations,  and  source  functions.  The 
time  domain  response  is  then  obtained  using  a  suitable  inverse  Fourier 
transform. 


2.2  Equivalent  Excitation 


It  is  desired  to  develop  an  equivalent  model  of  the  transmission 
line  problem  of  figure  1  .  Characteristics  of  the  model  should  include  a 
simplified  source  characterization,  ease  of  implementation  in  circuit  analysis 
codes,  applicability  to  nonlinear  loading,  and  a  relatively  easy  physical 
implementation  for  testing  purposes.  Equivalent  source  theory  will  be 
reviewed  briefly  before  an  alternate  coaxial  cable  coupling  model  is  devel¬ 
oped  . 


2.2.1  Thevenin's  and  Norton's  Theorems 


Both  Thevenin's  and  Norton's  theorems  deal  with  the  development  of 
equivalent  source  representations  of  a  linear  network,  N,  that  excites  an 
arbitrary  load  network,  N,  as  shown  in  figure  2.  Thevenin's  theorem  states 
that  the  network  N  may  be  replaced  by  an  equivalent  generator  consisting  of  a 
voltage  source  and  a  series  internal  impedance  such  that  the  terminal  response 
(current  and  voltage  at  T1,T2)  will  be  identical  to  the  original  network  (N 
and  N').  An  equivalent  current  source  representation  may  be  developed  through 
the  application  of  Norton's  theorem. 


Peskin  [5]  develops  Thevenin's  and  Norton's  theorems  subject  to  the 
following  assumptions: 

(a)  The  network  N  in  figure  2  is  linear. 

(b)  The  network  N  is  stable  over  the  full  range  cf  load  values 
possible . 

(c)  The  terminal  response  is  Laplace  transformable;  this  does  not 
require  N'  to  be  linear;  only  its  load  voltage  and  current 
must  obey 

oo 

J  I  f(t)  |  e°ldt<«>  , 

O 

for  some  real  number,  o. 
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Using  these  assumptions  and  the  node  relationships,  it  is  a  rela¬ 
tively  straightforward  matter  to  prove  that  the  network  N  may  be  replaced  by  a 
voltage  source  equal  to  the  open  circuit  voltage  at  terminals  T1 ,T2  and  a 
series  impedance  equal  to  the  input  impedance  of  N  looking  into  T1,T2.  The 
input  impedance  is  found  by  shorting  out  all  independent  voltage  sources  and 
opening  all  independent  current  sources.  The  Norton  equivalent  uses  a  current 
generator  equal  to  the  short  circuit  current  at  T1,T2  shunted  by  the  input 
admittance  (the  reciprocal  of  the  input  impedance). 


This  two-terminal  or  single  port  theorem  can  be  extended  to  multi- 
port  networks  if  we  require  m  independent  sources,  one  for  each  of  m  ports, 
and  if  we  define  an  impedance  (admittance)  matrix  composed  of  the  self  imped¬ 
ance  (admittance)  and  transfer  impedance  (admittance)  from  each  of  the  other 
ports.  The  source  voltages  for  a  Thevenin  equivalent  network  are  defined  by 


V?C  =  V. ,  I.  =  0,  j  =  I,  m 


(2) 


where  Vk  is  the  potential  at  port  k  and  I,  is  the  current  through  port  j. 
Similarly,  the  Norton  current  sources  are  tne  short  circuit  currents  at  each 
port  when  all  the  ports  are  shorted.  The  impedance  matrix  for  the  Thevenin 
representation  is  defined  by  the  equation 


z.j  =  Vi^j,Ik  =  0,VK^j 


(3) 


under  the  condition  that  all  internal  independent  sources  are  properly  shut 
off.  The  admittance  matrix  is  the  inverse  of  the  impedance  matrix.  The 
multiport  Thevenin  theorem  will  be  applied  in  the  next  section  to  develop  an 
equivalent  representation  of  the  coaxially  shielded  cable. 

2.2.2  Thevenin  Equivalent  Representation 

Modeling  the  coaxial  cable  response  with  a  Thevenin  equivalent  gen¬ 
erator  requires  computation  of  the  open  circuit  voltage  at  each  end.  Using 
equation  (1)  and  letting  Zf,Z2  approach  infinity  yields  the  following  open 
circuit  voltages: 

vo.  (0)  -  -  r v.  (p — - q  d%  -  k  j  i,  ® cosh  rq ' 11  d % 

i  sinh  n  {  sinhn 
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Figure  3  shows  the  equivalent  network.  The  impedance  matrix,  Z,  is 
[ K/tanh  H  -K/sinh  nl  (5) 


K/tanh  H 
Z  = 

K/sinh  n 

and  the  terminal  response  is  given  by 


-K/sinh  n 
-K/tanh  n 


_  v>  '■  v; 

v  rJ  N+  vf  =ZI+V> 


Using  the  impedance  relationship  for  the  load  voltage  and  current,  I1 
and  12  =  ^2^2'  ^e  l°ac*  voltage  is  found  from  equation  (6): 


V  =  Z 


-v,/z, 


w.  v+v> 


Equation  (7)  may  be  solved  for  V  as  follows: 

JmT  “If, 
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sinh  n  + 


K  cosh  nj 


sinh  n  + 


K  cosh  n  j 


f 

V 

1 

l 


Figure  3.  Equivalent  network  for  coaxial  cable. 
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where  u  is  the  identity  matrix.  After  substitution  of  equation  (4)  for  the 
source  term,  Vs,  manipulation  of  equation  (8)  yields  equation  (1),  demonstra¬ 
ting  the  correctness  of  the  equivalent  network. 

The  advantages  of  the  equivalent  representation  are  that  the  dis¬ 
tributed  source  is  replaced  by  two  point  sources  and  the  excitation  is  sep¬ 
arated  from  the  transmission  line.  This  means  that  standard  transmission-line 
models  available  in  most  circuit  analysis  codes  may  be  used  in  conjunction 
with  two  voltage  sources  to  model  the  coaxial  cable.  Similarly,  from  a  simu¬ 
lation  standpoint,  an  injection  system  may  be  developed  using  two  independent 
voltage  sources,  one  for  each  end  of  the  cable.  Therefore,  the  objective  of 
developing  a  simplified  representation  for  the  EMP  coupling  to  a  coaxial  cable 
has  been  achieved.  The  one  area  not  fully  explored  is  the  correctness  of  the 
equivalent  representation  when  nonlinear  loads  appear  at  the  terminals  in 
place  of  Z1  and  Z2.  Nonlinear  loading  is  studied  in  the  next  section. 

2.2.3  Nonlinear  Terminations 


It  was  demonstrated  that  the  equivalent  Thevenin  representation  for 
the  coaxial  cable  produced  results  identical  to  the  distributed  excitation 
when  linear  terminations  are  used.  While  the  single  port  representation  is 
valid  for  nonlinear  loads,  nothing  has  been  done  to  validate  the  multi-port 
representation  of  distributed  networks  with  nonlinear  loads.  In  the  following 
analysis  it  is  assumed  that  the  transmission  line  is  linear  and  that  the 
potentials  created  because  of  the  EMP  excitation  do  not  exceed  the  dielectric 
breakdown  voltages  for  the  cable.  Note  that  this  does  not  preclude  the  inclu¬ 
sion  of  nonlinear  effects  on  the  exterior  of  the  cable,-  only  its  internal  line 
must  be  linear.  If  internal  breakdown  is  possible,  then  the  model  would  have 
to  be  further  partitioned  to  isolate  the  nonlinear  regions. 

Obtaining  a  closed  form  solution  for  the  response  of  a  transmission 
line  with  nonlinear  loads  is  in  all  probability  impossible.  However,  if  a 
slight  restriction  is  placed  on  the  nature  of  the  terminations,  the  problem 
becomes  fairly  straightforward ,  and  useful  insight  may  be  gained  into  the 
actual  line  response.  As  illustrated  in  figure  JJ,  it  is  possible  to  very 
closely  approximate  a  nonlinear  load  characteristic  with  a  set  of  linear 
curves.  The  piecewise-linear  approach  will  be  used  to  study  the  validity  of 
the  equivalent  transmission  line  representation  under  nonlinear  loading. 
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Figure  Piecewise-linear  representation  of  nonlinear  load. 

Solving  for  the  response  with  nonlinear  loads  simply  involves  a  series  of 
linear  solutions,  with  the  loading  over  a  specific  time  period  being  deter¬ 
mined  by  the  load  voltage  or  current.  The  only  complication  in  the  solution 
is  the  inclusion  of  initial  condition  effects  at  transition  points. 


An  example  of  how  a  piecewise-linear  solution  is  used  is  presented 
in  figure  5.  Here  the  response  of  a  circuit  loaded  by  the  simple  square  law 
device  of  figure  is  computed.  As  the  load  current,  I,  reaches  a  transition 
value  of  the  piecewise-linear  model,  a  new  linear  equation  must  be  solved.  In 
practice,  only  a  finite  number  of  segments  will  be  used  to  represent  the  non¬ 
linear  load.  However,  if  computation  load  is  not  a  consideration,  then  a  very 
large  set  of  linear  curves  may  be  used  to  achieve  any  desired  accuracy.  The 
current  response  shown  in  figure  5  was  calculated  with  only  two  linear  re¬ 
gions.  The  piecewise-linear  response  compared  within  one  percent  with  the 
response  computed  with  the  MICROCAP-II  circuit  analysis  code. 


The  solution  for  the  EMP  response  of  the  coaxial  cable  given  in 
equation  (1)  assumed  that  the  initial  conditions  were  zero.  The  effect  of  the 
nonzero  initial  conditions  at  different  points  in  time  is  seen  in  the  trans¬ 
forms  of  the  differential  equations  describing  the  transmission  line 
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Figure  5.  Piecewise-linear  solution. 


-sr  =  -(R+jo)L)I  +  LI(to)  +  V 


-r—  =  -  (G  +  jcoC)V  +  CV(t )  +  I 


The  effect  of  the  initial  conditions  is  to  add  constant  source 
terms,  LI(tQ)  and  CV(tD),  to  the  EMP  sources,  Vs  and  Ig.  Defining  two  new 
source  terms, 


and 


(10) 


V  =V  +  L I  (t  ) 

s  s  o 

I  =  I  +CV(t); 

s  s  v  o' 


the  transmission-line  response  is  given  by  equation  (1)  with  the  new  source 
terms.  Therefore,  from  a  computational  standpoint,  a  new  Thevenin  equivalent 
generator  is  needed  to  model  the  line  after  a  transition  has  occurred  in  one 
of  the  terminations .  The  impedance  matrix,  Z,  remains  unchanged,  but  the 
sources  at  each  end  are  now 


Voc  (°>  =  •  j  Vs  © sin-  r(-  ~  ^  dS  -  K  j  l'  £)  ??shrG:.P 

J  sinh  n  J  sinhn 

o  o 
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Therefore,  the  criterion  for  determining  the  validity  of  the 
Thevenin  equivalent  network  with  nonlinear  loading  is  whether  or  not  it  pro¬ 
duces  the  open  circuit  voltages  defined  in  equation  (11)  after  load  transi¬ 
tions.  The  question  arises  not  because  of  the  distributed  EMP  sources,  Vg  and 
Ig,  but  rather  because  of  the  initial  conditions,  I(tQ)  and  V(tQ).  Obviously, 
the  transmission-line  current  and  voltage  distributions  are  different  for  the 
actual  EMP  excitation  and  its  Thevenin  equivalent.  If  the  proper  open  circuit 
voltage  is  produced  after  time  tD,  then  the  Thevenin  equivalent  network  i3 
valid  even  when  nonlinear  loads  are  used. 

One  method  of  verifying  the  correctness  of  the  open  circuit  volt¬ 
ages  would  be  to  (1)  compute  the  complete  line  response  under  both  types  of 
excitation,  (2)  transform  the  response  into  the  time  domain  to  obtain  the  ini¬ 
tial  condition  distribution  along  the  line,  and  (3)  solve  for  the  open  circuit 
voltages  with  the  distributed  initial  conditions,  V(tQ)  and  I(t0),  and  EMP 
source  terms,  Vg  and  Ig.  Unfortunately,  it  is  very  difficult  to  verify  open 
circuit  voltages,  even  for  a  simple  circuit.  A  similar  approach  was  used  to 
analyze  a  two-section  lumped  parameter  transmission-line  model  [6].  Based  on 
the  complexity  encountered  with  only  a  simple  circuit,  it  is  not  likely  that  a 
general  solution  can  be  obtained.  Therefore,  a  different  approach  is  needed. 

Assume  that  the  terminal  responses  of  a  coaxial  cable  are  as  shown 
in  figure  6.  Let  there  be  a  nonlinearity  in  the  termination  at  x  «  o  such 
that  the  impedance,  Z1 ,  changes  due  to  the  voltage  at  time  tQ.  Next,  assume 
that  the  open  circuit  voltages  from  time  tQ  on  are  as  shown  in  figure  7.  Note 
that  these  open  circuit  voltages  are  different  from  the  voltages  starting  at 
time  tQ  due  to  the  loading  from  time  zero  to  tQ.  These  open  circuit  voltages 
may  now  be  used  to  compute  the  response  of  the  line  after  the  nonlinear  tran¬ 
sition  at  time  tQ  as  shown  in  figure  6  by  the  dashed  waveforms. 
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Using  contradiction,  we  will  now  prove  the  validity  of  using  a 
Thevenin  equivalent  network  to  model  the  response  of  a  transmission  line  with 
nonlinear  loading.  In  section  2.2.2  it  was  shown  that  the  response  of  the 
Thevenin  network  would  duplicate  the  line  response  with  linear  loads  (the 
solid  curves  in  fig.  6).  Now  assume  that  the  open  circuit  voltages  of  the 
equivalent  circuit  after  time  tQ  are  not  the  same  as  figure  7.  This  is  the 
only  way  the  equivalent  network  will  not  provide  the  correct  response  after 


the  transition  at  time  t. 


We  now  have  two  different  networks  for  times 


greater  than  tQ,  one  with  the  correct  open  circuit  voltages  from  the  actual 
transmission  line  and  one  with  incorrect  open  circuit  voltages  from  the  origi¬ 
nal  Thevenin  equivalent  network.  The  impedance  networks  are  identical  because 
of  the  assumption  that  the  line  is  linear  and  operating  below  breakdown.  In 
section  2.2.2  it  was  shown  that  a  Thevenin  equivalent  network  will  provide 
correct  terminal  responses  for  any  linear  loading.  Now  reconnect  the  original 
terminations,  Z1  and  Z2,  to  the  new  equivalent  network.  Since  it  was  assumed 
that  the  open  circuit  voltages  from  the  original  Thevenin  model  are  incorrect, 
the  response  with  the  original  loading  is  also  incorrect.  This  violates  the 
work  of  section  2.2.2  that  proved  that  the  Thevenin  equivalent  would  provide 
the  proper  response  for  all  time  with  a  linear  load.  Therefore,  the  assump¬ 
tion  must  be  false  that  the  open  circuit  voltages  produced  by  the  equivalent 
network  are  different  from  the  open  circuit  voltages  produced  by  the  actual 
distributed  excitation. 

An  alternate  approach  for  validating  the  Thevenin  equivalent  net¬ 
work  can  be  developed  that  uses  the  uniqueness  of  the  solution  for  the  trans¬ 
mission  line  differential  equations. 

Since  the  differential  equations  governing  the  transmission  line 
are  linear  and  time  invariant,  the  solution  for  a  given  input  and  initial 
state  is  unique.  This  implies  that  the  input  required  to  produce  a  specified 
output  from  an  initial  condition  set  is  also  unique. 

The  open  circuit  voltages,  Vg,  required  to  produce  a  given  terminal 
response,  and  V2(tQ),  starting  from  time,  tQ,  are  found  from  equation 
(8)  to  be 


0 

l/Zj 


VA> 

YA> 


(12) 


where  the  voltages  ( tQ )  are  the  Fourier  transforms  of  the  time  domain  wave¬ 
forms  from  time  tQ  on.  Since  the  actual  network  response  and  its  Thevenin 
equivalent  produce  identical  terminal  responses  with  linear  loads,  Z i  and  Z2, 
the  open  circuit  voltages  from  time  tQ  must  be  identical.  Otherwise,  more 
than  one  set  of  inputs  will  provide  the  same  outputs.  This  is  not  possible 
with  a  linear  line.  The  effective  open  circuit  voltages  starting  at  an  arbi¬ 
trary  time  tQ  due  to  the  impressed  sources  and  the  initial  conditions  must  be 
the  same  for  the  actual  line  with  distributed  excitation  and  its  Thevenin 
counterpart  with  only  two  source  terms.  Therefore,  both  lines  will  provide 
the  same  terminal  response  after  a  termination  transition  at  a  time,  tQ. 


Equation  (12)  suggests  a  method  for  solving  for  the  response  of  a 
coaxial  line  with  nonlinear  terminations  using  Fourier  transform  techniques. 
The  load  nonlinearity  is  assumed  to  be  piecewise-linear .  The  steps  to  solve 
for  the  nonlinear  response  are  as  follows: 

(a)  Compute  the  complete  terminal  response  with  the  nominal  values 
for  the  terminations.  This  is  done  in  the  frequency  domain 
and  transformed  into  the  time  domain. 

(b)  Determine  the  time,  t0,  where  a  transition  will  occur  in  a 
termination. 

(c)  Compute  the  effective  open  circuit  voltage  with  a  start  time 
of  tQ.  These  voltages  are  found  by  computing  the  Fourier 
transforms  of  the  terminal  responses  multiplied  by  a  unit  step 
function  starting  at  tQ.  Alternatively  this  may  be  accom¬ 
plished  by  convolution  of  the  frequency  response  at  the  termi¬ 
nals  with  a  unit  step  defined  by 


(13) 


(d)  Solve  for  the  terminal  response  after  the  transition  (t  2  tQ) 
using  equation  (8),  the  source  voltages  from  step  3,  and  the 
new  terminations. 

This  is  a  simple  approach  to  solving  a  complex  problem.  This 
technique  should  be  computationally  superior  to  finite  difference  techniques 
when  the  line  lengths  are  long  and  many  cells  would  be  needed  to  model  the 
transmission  line. 


3.  RESPONSE  OF  A  MULTICONDUCTOR  TRANSMISSION  LINE 

The  equivalent  source  representation  developed  in  section  2  for  simple 
coaxial  cables  will  now  be  extended  to  multi conductor  transmission  lines. 
After  a  brief  review  of  raulticonductor  transmission  line  theory  in  section 
3.1,  the  equivalent  source  representation  is  developed  in  section  3*2.  Sever¬ 
al  implementation  techniques  are  reviewed  in  section  3.3. 

3.1  Distributed  Excitation 

Consider  the  lossless  multiconductor  transmission  line  shown  in 
figure  8.  The  line  consists  of  N  +  1  conductors,  with  conductor  N  +  1  being 
the  reference  conductor.  The  length  of  the  line  is  1,  and  its  electrical 


Figure  8.  Multiconductor  transmission  line. 

characteristics  are  determined  by  the  capacitive  coefficient  matrix  C  and  the 
inductive  coefficient  matrix  L.  The  line  is  described  by  the  following  matrix 
partial  differential  equations  [7]: 

3v(x,t)  _  3i(x,t) 

3x  3t 


3i(x,t)  _  ^  3v(x,t) 
dx  dt 


when  no  source  terms  are  present.  Using  the  Laplace  transform  these  equations 
become : 


['I- 


where  the  dependence  of  7  and  I  on  x  is  implicit;  V  and  I  are  n  *  1  vec¬ 
tors;  C  and  £  are  n  x  n  matrices  and  5  is  a  zero  matrix  of  order  n  x  n. 

Differentiating  equation  (14)  and  substituting  yield  a  set  of 
decoupled  second-order  differential  equations: 

rvi  .  rtc  o  l  fvi 


Ct  111 


It  is  possible  for  lossless  multiconductor  cables  to  define  a  transformation 
matrix  T  such  that 


s2?'1  CLT  =  y2  , 


(17) 


Figure  9.  Generalized  point  source  for  multiconductor  transmission  line. 


where  Y  is  an  n  x  n  diagonal  matrix  with  elements 


r*.-j 

0,i  *j 


(18) 


and  Y-|  an  eigenvalue  of  equation  ( 1  6 ) . 


Then  using  the  transformations  V  =  TV  and  I  =  TI  ,  equation  (16) 
a  set  of  2N  decoupled  differential  equations  with  a  solution  given  by 


becomes 


(19) 


where 

-  -1  -k-1  -  “  --1 

Zc  =  s  E  TyT 


(20) 


is  the  characteristic  impedance  matrix,  U  is  the  identity  matrix,  and  a+ 
and  a~  are  n  x  1  vectors  determined  by  the  terminating  impedances  and  the 
source  terms  exciting  the  line.  Tesche  et  al  [7]  provide  a  solution 
for  a+  and  i-  when  the  line  is  terminated  in  impedance  matrices  at  x  *  0 
and  x  =  1  ,  respectively,  and  in  a  generalized  point  source  at  xg,  as  shown  in 
figure  9.  Of  interest  here  is  the  terminal  response  of  the  line.  The  voltage 
and  current  at  the  ends  are  related  by  the  terminating  matrices  as  follows: 


V(0) 

V(l) 


and  from  Frankel  [8], 


(21) 


(22) 


where 

n=[VZ]  [^i'^c]  (23) 

is  the  rjeflection  coefficient  associated  with  load  impedance  Zi , 
and  6-  and  6  are  source-related  terms.  The  operator  in  the  above 

equations  arises  from  the  fact  that  the  elements  of  the  matrices  are  also 
matrices.  The  matrices  are  first  multiplied  as  if  they  contained  normal 
scalar  elements.  Then  the  resulting  matrix  operations  for  each  element  are 
carried  out. 


The  source  terms  5  and  5+  have  been  presented  [9]  for  a  general¬ 
ized  point  source  (fig.  9)  as 


(24) 


where  xg  is  the  location  of  the  source.  If  a  distributed  source  excites  the 
line  then  equation  (22)  must  be  integrated  over  the  source: 


1 

5*  =  }st©d5  , 

0 

and  the  voltage  and  current  sources  (Vg  and  Ig)  are  per 
quantities . 


(25) 


unit  length 


3-2  Equivalent  Excitation 

As  with  the  single  coaxial  cable,  it  is  desired  to  develop  an 
equivalent  representation  that  readily  lends  itself  to  circuit  analysis  and 
point  source  simulation.  The  Thevenin  equivalent  network  for  the  multiconduc¬ 
tor  cable  is  a  straightforward  extension  of  the  coaxial  cable  development. 
The  source  at  each  end  must  be  an  n-dimensional  vector  to  satisfy  the  most 
general  conditions.  The  2n  port  impedance  matrix  now  requires  2n  *  2n  terms 
to  account  for  the  self  and  transfer  impedances  of  all  the  ports.  The  cor¬ 
rectness  of  the  multiconductor  Thevenin  equivalent  network  will  be  examined 
first.  Then  the  impact  of  a  restriction  on  the  Thevenin  source  terms  will  be 
investigated.  The  purpose  of  the  source  restriction  is  to  further  simplify 
the  model . 


The  response  of  the  shielded  multi  conductor  cable  to  the  Thevenin 
equivalent  source  excitation  may  be  obtained  from  equations  (21)  through 
(25).  The  layout  of  the  Thevenin  equivalent  network  is  given  in  figure  10, 
where  Vg(0)  and  Vg(1)  are  the  wire  to  reference  open-circuit  voltage  vectors 
at  x  *  0,  x  «  1,  aue  to  the  distributed  EMP  excitation: 


3.2.1  Thevenin  Source  Terms 


(26) 


The  open  circuit  voltages  Voc(o)  and  Voc(1)  are  found  from  equa¬ 
tions  (21)  through  (23)  and  equation  (26)  if  we  let  the  terminating  impedances 
between  each  wire  and  the  reference  conductor,  as  well  as  the  wire-to-wire 
impedances,  approach  infinity  (open  circuit).  The  impedance  matrix  Zi  may  be 
found  from  the  definition 
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Vv/H=0Vk*j 


(27) 


yielding 


Z  = 

*j 


r  i=j 
1  0,i*j  • 


(28) 


“  IT  T 

The  source  vector  [6  ,6  ]  in  equation  (26)  is  independent  of  the  terminating 
matrices. 


The  elements  of  equations  (21)  through  (23)  that  depend  on  the 
termination  matrices  are  of  the  form 


2,  [  U  -  ?,]  =  22|  [Z,  -  2  ]  Ze 


(29) 


and 

r-[Vzj'  [V2J  .  (30) 

In  the  limit  ( Z ^  defined  by  eq  (28)),  it  is  found  that  equation  (29)  ap¬ 
proaches  2Z  .  The  reflection  matrix,  ,  approaches  the  unit  matrix,  0,  under 
open  circuit  conditions. 


Figure  10.  Equivalent  source  representation  for  multiconductor  transmission 


Using  equations  (21 )  and  (22)  to  define  the  terminal  voltages  and 
taking  the  limit  as  the  terminating  impedances  approach  infinity  yields  the 
following  open  circuit  voltages  at  the  terminals: 


— oc 

V  (o) 

—  OC 

V  (l) 


Te'7*  T'1 


Te"7*  T  1 


with  the  source  terms  and  5^  as  given  in  equation  (26). 

3-2.2  Terminal  Response 

The  terminal  response  for  the  Thevenin  equivalent  network  is  found 
by  superposition  of  the  response  from  the  source  at  each  end  of  the  cable. 
The  source  voltages  given  in  equation  (31)  are  used  in  conjunction  with  equa¬ 
tion  (24)  to  define  the  source  terms  fil  and  6*  for  the  current  responses  given 
in  equation  (22).  The  current  source,  Ig,  in  equation  (24)  is  zero  because 
only  a  voltage  source  vector  is  used  at  each  end  of  the  cable.  Summing  the 
source  terms  due  to  the  voltage  sources  at  xg  =  0,1  yields 
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\l  g+1 

where  L  d  d  J  is  the  distributed  source  term  vector.  This  shows  that  the 
Thevenin  equivalent  network  provides  the  same  excitation  at  the  terminals  as 
the  distributed  excitation.  The  terminal  response  is  given  by  equations  (21) 
and  (22). 
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3.2.3  Nonlinear  Terminations 


The  nonlinear  termination  study  for  a  single  coaxial  line  (sect. 
2.2.3)  is  directly  applicable  to  multiconductor  lines.  The  only  real  question 
was  whether  or  not  Thevenin  equivalent  circuits  of  transmission  lines  are 
valid  with  nonlinear  terminations.  A  multiconductor  cable  simply  requires  a 
more  complex  source  and  impedance  representation.  Therefore,  the  work  of 
section  2.2.3  implies  that  equivalent  circuit  representations  are  valid  for 
multiconductor  transmission  lines  with  nonlinear  terminations. 

3.2.4  Coupling  Distribution 


In  the  previous  sections  it  wts  shown  that  an  equivalent  source 
network  of  2N  independent  voltage  terms  will  provide  the  same  terminal  re¬ 
sponse  as  the  distributed  excitation.  The  distribution  of  the  shield  coupling 
parameters  among  the  N  conductions  will  now  be  examined  in  an  attempt  to 
further  reduce  the  source  requirements. 

The  distribution  of  the  coupling  through  the  shield  will  determine 
the  source  terms  Vg  and  Ig  in  equation  (22).  Unfortunately,  little  work 
exists  on  the  nature  of  the  coupling  distributions.  As  pointed  out  in  refei — 
ence  [8],  most  shield  coupling  studies,  both  experimental  and  theoretical, 
consider  coaxially  shielded  cables  with  only  one  conductor.  Frankel  [8] 
suggests  measurement  techniques,  but  does  not  provide  any  insight  into  the 
distribution. 

In  Tigner  et  al  [9],  a  theoretical  treatment  is  proposed  for  the 
distributions  of  the  capacitive  and  inductive  coupling  parameters;  the  common 
mode  inductive  and  capacitive  sources  are  distributed  over  the  conductors 
using  coupling  ratio  vectors  c  and  1  3uch  that 


£c=i,  £i=i. 

i=l  i=l 


The  requirement  of  equation  (33)  appears  plausible  for  the  capaci¬ 
tive  coupling  parameter  since  in  effect  the  measured  common  mode  current 
source  is  really  a  set  of  N  current  sources  in  parallel,  each  driving  one  of 
the  wires.  However,  the  requirement  that  the  N  multi conductor  voltage  sources 
add  up  to  the  common  mode  voltage  sources  appears  to  be  inconsistent  since  the 
N  voltage  sources  are  also  in  parallel,  not  in  series  as  equation  (33)  sug¬ 
gests.  Thi3  appears  on  average  to  divide  the  voltage  source  amplitude  by  the 
number  of  conductors.  Evidence  of  this  improper  scaling  is  seen  in  the 
results  presented  by  Paul  [10]  for  a  comparison  of  the  theoretical  model  with 
experimental  data.  A  very  large  shield  transfer  impedance  of  1  ohm/meter  was 
derived  empirically  for  a  37-conductor  cable  with  a  tape-wrapped  shield.  It 
is  felt  that  this  large  shield  transfer  impedance  is  an  artifact  of  the 
requirement  of  equation  (33). 


Unfortunately,  no  other  models  or  experimental  data  have  been  dis¬ 
covered  in  the  literature.  The  following  model  development  attempts  to 
overcome  the  difficulties  discussed  above.  Only  the  inductive  source  is 
considered  in  this  study.  Experimental  data  suggest  that  this  is  the  dominant 
mode  of  coupling  for  most  good  quality  shields.  Figure  11  presents  a  section 
of  multiconductor  line  with  an  overall  shield.  The  transfer  impedance  of  the 
shield  is  Zab  and  the  external  shield  current  is  Ig  such  that  the  potential 
along  the  inside  of  the  shield  is 

Vs  =  ZabIsAx  •  (3*0 

It  is  assumed  that  the  external  current,  Ig,  is  evenly  distributed 
around  the  cable.  The  line  integral  of  the  electric  field  is  related  to  the 
surface  integral  of  the  magnetic  field  by  Stokes  theorem: 


J*  E  •  dl  =  -jco  j*  B  •  dA  , 


where  A  is  the  area  enclosed  by  the  line  integral.  Assuming  that  the  dis¬ 
tances,  Ax,  and  the  wire-to-shield  separation  are  small,  the  line  integral  of 
the  electric  field  may  be  replaced  by  the  potentials  shown  in  figure  11: 


V.  (x)  -  V.  (x  +  Ax)  +  Vs  =  -jcojB«dA  * 


where  the  field  along  the  conductor  is  zero  due  to  the  assumption  of  perfectly 
conducting  wires.  Substituting  equation  (3*0  for  Vg,  dividing  by  x,  and  tak¬ 
ing  the  limit  as  x  approaches  zero  yields 


dV(x) 

dx 


=  -j(D  J 


B  •  dl  -Z  ,1 

ab  s 


with  h  being  the  separation  distance  of  the  wire  from  the  shield.  The  mag¬ 
netic  field,  B,  is  the  total  field,  incident  and  scattered.  The  scattered 
magnetic  field  may  be  related  to  the  wire  currents  through  self  and  mutual 
inductances.  The  incident  magnetic  field  is  the  field  present  if  the  internal 
wires  are  not  present.  However,  for  frequencies  of  interest  (£100  MHz)  and 
typical  cable  dimensions  (radius  <2  cm)  the  cable  without  wires  acts  as  a 
waveguide  well  below  the  cutoff  frequency.  Therefore,  the  incident  magnetic 
field  is  zero,  yielding  the  result  that  the  voltage  source  on  each  wire  is 
identical  and  equal  to  the  product  of  the  transfer  impedance  and  the  external 
current.  This  does  not  mean  that  the  open  circuit  voltages  at  the  ends  of  the 
cable  are  equal,  simply  that  the  per  unit  excitations  are  equal. 
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Figure  11.  Multiconductor  transmission-line  segment. 


Some  insight  may  be  obtained  into  the  open  circuit  voltage  response 
if  a  slightly  different  solution  representation  is  used  instead  of  equations 
(21),  (22),  and  (2*0 .  The  chain  matrix  formulation  provided  by  Paul  [10]  is 
as  follows: 
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and 


(40) 


Ol2(l)  =  -ZcTsinh  (jcoryi)  T 1 


<D22(1)  =  T  cosh  (jcoyl)  T  . 


The  open  circuit  voltage  may  be  found  from  equations  (36)  through 
(38)  if  we  simply  let  I(o)  and  1(1)  equal  zero. 


Now  consider  the  simple  condition  where  the  current  source  Ig  is 
negligible,  and  also  assume  that  the  external  shield  current  is  constant  over 
the  length  of  the  cable.  Then  the  open  circuit  voltage  at  the  x  =  o  end 
becomes 
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The  open  circuit  voltage  at  the  other  end  is  identical  except  for 
the  sign.  If  the  dielectric  is  homogeneous  then  the  transformation  matrix  T 
is  the  unit  matrix  and  the  eigenvalue  matrix,  Y,  becomes 


l 


(42) 


where  v  is  the  propagation  velocity  of  the  line.  Then  the  open  circuit  volt¬ 
age  becomes 

V 

V(0)=- 

jco  sinh  j  —  (43) 

v 

which  means  that  the  open  circuit  voltage  is  constant  for  all  the  wires  if  the 
source  voltage,  \?g,  is  constant.  The  open  circuit  voltage  at  x  -  1  is  iden¬ 
tical  except  for  the  sign.  A  similar  result  may  be  obtained  for  the  current- 
source-only  condition,  except  the  sign  does  not  change  from  end  to  end. 


If  the  cable  is  not  homogeneous  then  there  may  be  distinct  eigen¬ 
values.  Since  the  eigenvalue  matrix  Y  is  diagonal,  equation  (41)  may  be 
rewritten  as 


and  the  open  circuit  voltage  vector  will  not  be  constant  even  if  Vg  is  con¬ 
stant.  The  amplitude  of  each  mode  is  proportional  to  its  propagation  veloc¬ 
ity,  ,  and  the  wire  potentials  are  weighted  sums  of  these  mode  voltages  with 
the  weights  determined  by  the  T  matrix  and  the  characteristic  impedance 
matrix,  Zc< 


Therefore,  under  the  assumptions  of  inductive  coupling  dominance 
and  a  homogeneous  dielectric  the  open  circuit  voltage  from  each  wire  to  the 
shield  will  be  equal.  This  means  that  a  single  source  may  be  used  in  series 
with  the  cable  shield  at  each  end  of  the  cable.  This  allows  a  significant 
reduction  in  the  complexity  of  a  Thevenin  equivalent  network  for  multicon¬ 
ductor  cables.  Now  only  two  sources  are  needed,  not  2N.  For  nonhomogeneous 
dielectrics,  the  variation  in  wire  to  shield  potentials  may  be  small  enough  to 
allow  use  of  only  one  source  at  each  end.  The  suitability  of  this  approxi¬ 
mation  will  depend  on  the  specific  characteristics  of  the  cable  under  study. 
Implementation  details  are  discussed  in  the  following  section. 

3-3  Equivalent  Excitation  Implementation 


There  are  at  least  three  implementation  techniques  that  can  be  used 
in  the  simulation  of  the  EMP  response  of  a  shielded  cable.  The  purpose  of 
discussing  implementation  techniques  here  is  to  provide  theoretical  insight 
into  each  method  and  to  discuss  the  advantages  and  disadvantages  of  each.  The 
discussion  is  limited  to  implementation  on  physical  systems  since  circuit 
analysis  modeling  is  relatively  straightforward. 

If  it  is  known  that  the  open  circuit  voltages  are  widely  different 
for  each  of  the  N  conductors,  then  the  only  implementation  possible  is  the 
development  of  independent  sources  for  each  of  the  2N  ports.  These  sources 
must  be  reasonable  approximations  to  ideal  voltage  sources.  This  requires  the 
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source  impedances  to  be  much  smaller  than  the  characteristic  impedances  of  the 
cable.  The  sources  are  placed  in  series  with  each  conductor. 

Of  most  interest  here  is  the  condition  discussed  in  section  2.2.4, 
where  the  open  circuit  voltages  are  all  approximately  equal  at  each  end  of  the 
cable.  Then  only  two  sources  are  needed,  one  for  each  end.  These  sources 
must  be  inserted  into  the  reference  conductor  (the  cable  shield).  The 
simplest  way  to  implement  these  sources  is  to  insert  very  low  resistance  in 
series  with  the  cable  shield  and  create  the  required  voltage  waveform  across 
these  resistors;  see  figure  12.  The  sources  must  provide  currents  defined  by 

V*)  =  voc(x)/R*  x  -  0.  1  .  (45) 

where  the  resistance,  R,  must  be  much  less  than  the  bundle-to-shield  impedance 
of  the  cable.  The  key  advantages  of  this  method  are  low  power  requirements 
and  ease  of  specifying  the  source  waveforms.  The  major  disadvantage  is  that 
the  shield  continuity  is  disrupted  by  the  insertion  of  the  resistors. 

A  tradeoff  may  be  made  between  the  advantages  and  disadvantages  of 
the  simple  resistive  insertion  technique.  The  integrity  of  the  shield  may  be 
preserved  if  a  segment  of  the  shield  is  used  instead  of  resistors  to  obtain 
the  low  source  impedance.  The  source  current  waveform  is  then  defined  by 

Ig(x)  =  VQC(x)/(Rs  +  jwLs)dx,  x  •=  0,  1  ,  (46) 

where  Rg  and  Lg  are  the  per  unit  length  shield  resistance  and  inductance,  and 
dx  is  the  length  of  the  source  segment.  Equation  (46)  assumes  that  the  shield 
excitation  is  dominated  by  the  resistive-inductive  transfer  impedance  and  that 
the  segment  length,  dx,  i3  short  compared  to  the  smallest  wavelength  of  the 
open  circuit  voltage.  The  preservation  of  the  shield  is  obtained  at  the 
expense  of  the  source  current  requirement.  Typical  shield  transfer  impedances 
range  from  1  to  10  milliohms  per  meter,  whereas  the  source  resistance,  R,  may 
be  as  large  as  1  ohm  in  some  instances  and  still  satisfy  the  requirement  of 
being  much  less  than  the  bundle-to-shield  impedance.  Therefore,  the  source 
current  requirements  may  go  up  as  much  as  an  order  of  magnitude  or  more.  The 
sources  can  be  connected  to  the  shield  segments  if  some  of  the  insulation  at 
the  ends  of  the  cable  is  removed  or  if  a  short  extension  is  added  to  the 
cable.  The  main  advantage  of  this  method  is  the  preservation  of  the  external 
shield  integrity.  The  primary  disadvantage  is  the  high  source  current 
required,  with  a  secondary  problem  being  the  difficulty  associated  with  the 
physical  connection  to  the  shield. 

The  third  technique  attempts  to  reduce  the  disadvantages  of  the 
second  technique  through  excitation  of  the  entire  cable  shield  as  shown  in 
figure  13.  This  technique  has  a  much  more  complicated  source  specification 
because  of  the  very  distributed  nature  of  the  coupling  through  the  shield. 
The  source  current  is  given  in  the  frequency  domain  by 


where  Vd(x)  is  defined  as  the  open  circuit  bundle  voltage  due  to  an  impulse 
function  excitation  (constant  frequency  spectrum).  An  analytical  expression 
for  Vd(x)  may  be  obtained  from  equation  (4),  with  the  voltage  and  current 
sources  being  functions  of  the  external  response  to  the  impulse  functions. 
This  analysis  is  similar  in  complexity  to  the  normal  EMP  analysis  of  the 
cable.  It  is  also  possible  to  determine  Vd(x)  experimentally.  The  advantages 
offered  by  this  technique  are  a  slightly  better  source  power  requirement  over 
the  shield  segment  technique  and  complete  maintenance  of  the  shield  integrity. 
The  disadvantages  are  (1)  a  more  complex  source  specification  and  (2)  the 
possibility  that  the  source  current  may  not  be  a  realizable  function  due  to 
nulls  or  peaks  in  Vd(x). 
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Figure  12.  Equivalent  source  implementation  using  resistors  in  shield, 
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Figure  13-  Equivalent  source  implementation  over  entire  cable  shield. 


4.  CONCLUSION 


Efficient  equivalent  source  representations  were  developed  for  coaxially 
shielded  cables  and  for  multiconductor  cables.  The  correctness  of  these 
equivalent  source  representations  was  validated  for  linear  and  nonlinear 
terminations.  The  distribution  of  shield  coupling  parameters  was  studied  for 
shielded  multiconductor  cables.  The  results  of  this  study  indicate  that  an 
approximate  equivalent  source  excitation  may  be  used  for  multiconductor 
cables,  thereby  reducing  the  source  requirements  even  further.  The  approxi¬ 
mate  equivalent  source  excitation  requires  only  two  independent  sources  to 
represent  the  EMP  response  of  the  cable.  The  simplicity  of  the  equivalent 
source  excitation  means  that  numerical  models  will  run  quicker,  and  system- 
level  stimulation  is  possible  with  a  lower  power  level  and  less  complex 
simulators . 
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